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We study an array of cigar-like Bose atom condensates confined in a cylinder and examine the 
competition between the dipole-dipole and the short range interactions. The system is effectively 
reduced to a one- dimensional boson one with a contact and inverse square interactions. We call 
this system the Calogero-Sutherland-Lieb-Liniger gas. The universal properties of the ground state 
are analyzed by the renormalization group theory. By using the bosonization techniques to the 
excluson gas, we calculate the non-universal exponent depending on the microscopic parameters. 
This exponent may be experimentally measurable. 
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One-dimensional many-particle systems have been a 
central topics in physics for four decades. Many ex- 
actly soluble models were proposed after Bethe solved the 
one-dimensional Heisenberg model by his famous ansatz 
Q, 0. Among them two kinds of the exactly models 
are cross-sectional: The particles interacting with con- 
tact or with inverse square potentials 0, 0. The 
one-dimensional cold alkali metal Bose atom gas, which 
may condense at zero temperature @, is a typical ex- 
ample of the former. Here, we call it the Lieb-Linigcr 
(L-L) gas for Lieb and Liniger first gave the exact so- 
lution for this model. The strong interacting limit of 
the L-L model, known as the Tonks-Girardeau (T-G) 
gas 0, has also reached in the ultracold Bose atoms 
8J. The model with inverse square interaction, called 
the Calogero-Sutherland(C-S) models play a fundamen- 
tal role in the contemporary theoretical physics [!j . The 
C-S model may be a universal Hamiltonian of some 
weakly disordered metals [T^j. The models may also be 
used to describe the edge excitations in the fractional 
quantum Hall effect 0|. In a recent work we P r0 ~ 
posed a possible realization of the C-S gas for the cold 
atom with a dipole-dipole interaction, which is based on 
a major progress in the cold chromium atoms 52 Cr i.e., 
this cold atom gas is turned into a Bose condensate ^| . 
It may also be realizable in the cold atom in tight waveg- 
uides with anyonic exchange symmetry [hH ]. 

The models with both the short range and long range 
interactions mentioned above have been applied to the 
faceting transition of the miscut crystal surface^!, 
by renormalization group analysis [T^. The low en- 
ergy effective theory of these models may be character- 
ized by the non-ideal excluson gas 0] since the quasi- 
particle excitations in these models obey exclusion statis- 
tics [lJJ |2(J, |21| . It is a Luttinger liquid theory [23, Hj] 
whose non-universal exponent is given by both the short 
range interaction strength and the exclusion statistics pa- 
rameter. 

Both interactions are of the comparative magnitude in 
the cold chromium atom gas [l^. This offers an oppor- 
tunity to study the intergradation from one to another 
since both interactions are adjustable |2l|. In this 
Letter, we propose a realization of such a model for the 




FIG. 1: (Color on-line) The cigar-like condensates are con- 
fined in a cylinder, (a) The magnetic dipoles; (b) The electric 
dipoles. The insets show the angle relation in deriving eq. © 



cold atoms in a cylinder. Fig. ^is a sketch map of our set 
up. An array of cigar-like Bose atom condensates is con- 
fined in a cylinder which is made by, e.g., a ring-shaped 
optical lattice [2(j. The atoms interact with a contact 
potential accompanied a dipole-dipole potential [2~3 ■ An 
electromagnetic field is generated by a stable current. If 
the dipoles is magnetic, their orientations are shown in 
Fig. [Ja) while the orientations are as in Fig. db) if 
they are electric. The arc coordinate is x and the verti- 
cal direction is y. The two-dimensional wave function of 
the many atom system is factorized, i.e., the y-direction 
wave function is approximated by its condensate wave 
function and the system effectively is one-dimensional. 
The effective atom density of the system may be written 

as p(r) = po(y) X)i=i ^( x ~ x i) where po(y), approximated 
by a constant po, is the atom density in a single cigar- like 
condensate and N the number of the condensates. 

The three-dimensional interacting potential of the 
atom cloud is generally given by V = \ J s drdr' p(r)V (r— 
r')p(r') where S denotes the cylinder with the thick lo. 
For simplicity, we consider a dilute gas limit in which 
whatever an optical lattice exists, the continuous model 
works. The dilute condition is given by tqN/L << 1 
where L is the perimeter of the circle and tq is the az- 



2 




g 



FIG. 2: (Color on-line) The renormalization group flows for 
g=0.75 (orange), 0.2 (green), (red), -0.1 (blue), -0.25 (yel- 
low) and -0.6 (black). The arrows show the flow directions. 
The filled spots are the stable fixed points and the empty are 
the unstable fixed points. 



imuthal width of a condensate. The interaction Vfe 
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— for a is the s-wave scattering length. 



d is the dipolar moment of the atom. For the factorized 
wave function, the Hamiltonian may be reduced to a one- 
dimensional one 
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where G = 2d 2 pn, U\ oc h , a is the one-dimensional re- 
duction of U, whose detailed version depends on the con- 
straint potential. r c — > is the character size of the con- 
tact potential. S± = and the sign =p corresponds to 
the magnetic or electric dipoles in Fig^a) or (b). The 6± 
term in the Hamiltonian is not important in the thermo- 
dynamic limit because it is of order O(l) while the ground 
state energy of the system is an extended quantity. We 
call this system the L-L-C-S gas since it combines the C-S 
and L-L gases. In the following, we use the dimensionless 
coupling constants g = -p-G, u = U\, and v = u — g 
with t c /Iq fixed when r c and Iq — ► .In the system our 
are concerning, g > while u may be a real number. 

To understand the universal properties of the ground 
state, we employ the renormalization group theory. This 
has been built for models with both contact and long 
range interactions |T^ . For this inverse square model in 
the one-dimension, the renormalization group equation 
reads 
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v - v + g (3) by [| e x 



FIG. 3: (Color on-line) The fixed point lines. The 
solid(black) curve is the stable fixed points. In the red curves 
a pair of different g are corresponding to the same unstable 
fixed point u* . And in the green, the unstable fixed point for 
g > 0.75 shares the same u* with the stable fixed point. 



while the long range coupling constant g is a marginal op- 
erator. Here I = In N is the running parameter. The flow 
pattern 4r and the fixed point values v* are shown in the 
insets of Fig. |2]and Fig. 02 respectively 0,0. Since 
g > in the present model, one may write g = A(A — 1) 
for a real A. The fixed point value v* for a given g equals 
to A. In a pure C-S model, it is known that A serves as 
the exclusion statistics parameter 21] . This means that 
the fixed point line of this model is determined by exclu- 
sion statistics of the system. To understand the physics 
directly from the microscopic parameters, we redraw the 
flow and fixed point line in u— ii and g — u* in Fig. |2Jand 
Fig. 13 respectively. In Fig. [3 we see that there are two 
fixed points in a flow line for a given g > —0.25. The right 
hand side one is stable and the left one is unstable. Each 
fixed point value u* corresponds to two g± = u* ± y/u*. 
If < u* < 0.25, both fixed points are unstable. For 
u* > 0.25, one of g+ is unstable while g_'s is stable. Fig. 
13 clearly shows this fact. It is seen that there is no fixed 
point with u* < or g < -0.25. If g < -0.25, the at- 
traction between condensates (called particles hereafter) 
is so strong that the on-site repulsion can not withstand 
the collapse. This collapse, however, will not happen for 
the system we are considering, in which g > 0. 

If we only concentrate on the short range interaction, 
the above renormalization group analysis for the ground 
state properties may be complete. However, since the 
long range interaction may lead to consequences with 
the exception of the above renormalization group anal- 
ysis, we have to carefully examine the ground state. 
The renormalization group flows show that the pure C- 
S model with u = is driven to u — > — oo. We ex- 
amine the meaning of this running now. There are 
two real solutions for g = A(A — 1) if g > —0.25, i.e., 
A± = 0.5(1 ± \/I + 4g). The ground state energy is given 



Eli(2j-A r -1) 2 - For -0.25 < g < 0, 
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< A_ < 0.5 and the ground state wave function of the 
C-S model is given by 



c-s 



n 



n(xi - Xj) 



(4) 



This means that although the interaction is attractive, 
there is still a mutual repulsion between particles caused 
by the quantum zero motion such that the system is still 
stable against the collapse. Since the contact potential is 
less singular than the 1/x 2 potential, for —0.25 < g < 0, 
no matter what is the magnitude of u, the ground state 
wave function exactly given by Reflected in Fig. 

121 , there is a small region between the red curve and 
the yellow curve at u = (the black short stub), the 
system may be characterized by the C-S model. We have 
discussed such a situation in our previous work |T^ and 
do not concern it here because in our model g > 0. 

For g > 0.75, A_ < —0.5, and Q is not square inte- 
grable, the ground state wave function is given by [{| 




FIG. 4: (Color on-line) The effective Luttinger liquid param- 
eter A e //(the green surfaces) varies as A and Ui. The unit of 
U-i is hv s . (a) < A < 0.5; (b) A > 1.5. 
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with A+ > 1.5. Hence, the system will not run to u — 
— oo from the unstable fixed points but is characterized 
by the thick vertical line at u — in Fig. [21 

For < g < 0.75 (-0.5 < A_ < 0), there is a ground 
state (0J which is square integrable. The renormalization 
group flow indeed will drive the system to a state with 
u = — oo for a negative u fluctuation. Howe ver, at exact 
u = 0, because of the symmetry A <-> —A 28] lg, the 
system is equivalent to that with < A < 0.5. This 
character of the system is shown by the vertical dash line 
in Fig. 

According to the above analysis, we have an overall 
view to the ground states of the system. We see that 
there is no ground state with 0.5 < A < 1.5 in the present 
model. This is because the non-interaction particles are 
boson, the system has to back to the free boson when g = 
u = 0. If applying this model to fermions, A_ solutions 
have to be abandoned because A — ► 1 as g — > 0. 

In the cold atom context, many system parameters 
may be exactly controlled. Therefore, one should do mi- 
croscopic calculations to understand many non-universal 
properties. Here we use an effective theory early devel- 
oped in 01 1 t° study some low energy behaviors of the 
system. The theory was called the excluson gas theory. 
Due to the symmetry A <-* —A, we focus on A > 0. 

We first consider A ^ 0. The pseudo- Fermi- momentum 
is defined by kp = y/2mp for the chemical potential p 
. Its value is fixed by J_^ F dkp(k) — = do for No is 
the one-dimensional particle number in the ground state. 
The particle density p(k) = 1/(2ttXH) if \k\ < kF and if 
\k\ > kp. This gives kp = nXhdo and p = ^(irXhdo) 2 ■ 
The ground state energy can be written as, in the thermo- 

2 X 2 d~l 



dynamic limit, E\/L = ^ /_£ dkp(k)k 2 = 



the low-lying excitations of the system are: The density 
fluctuations with the sound velocity v s — vp — kp/m; 
Adding extra number M = N — No of particles to 
the ground state; And creating a persist current by the 
Galileo boost k — > k + irJh/L in the ground state. To 
obey the periodic condition, on requires (— 1) J = ( — 1) M 
for integers M and J. If we consider only these low-lying 
excitations, the non-ideal excluson gas may be described 
by a bosonized effective Hamiltonian 0] 
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where vn = Xv s and vr = 
tinger liquid relation |22l |. 



v s /X are the standard Lut- 
The boson operators b q and 



b q are related to the density fluctuating operators p q 



by b q = y/2-K/qL p q + \ b q = y/2%/qLp q )f with ± de- 
noting the right and left Fermi points. According to 
the periodic boundary condition, a consistent choice for 
M = Mr+M l and J = J R +J L is M R = J R , M L = -J L . 
Using the Bogoliubov transformation, the Hamiltonian 
can be easily diagonalized 



H B = ^2 v s q{a\a q 
q >o 



a\a q ) + ^[v N M 2 + vjJ% (7) 
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where the Bogoliubov boson is given by a q = b q cosh (po — 

The renormal- 

U 2 \^ 2 , v N = 
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ized velocities are v 
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to x~ Ae// so that the momentum distribution is given by 
n(p) ~ C\ — C*2Sgn(p — Pf)\p ~ pf^ X " ss ■ In the limit 
U\ — > oo, the system goes an insulator state because 
x~ x °ff — > as X e ff — > oo. 




FIG. 5: (Color on-line) The sketch map of the image of the 
time-of- flight. Circles represent the positions of the atoms 
with the fastest radial velocity at the different flight time. 
The cross points between each blue solid lines and the circle 
is the position of these atoms. The red dotted lines are the 
reference radial lines. 



v s e 2 ^°, and vj = v s e 2 ^ with the controlling param- 
eter ipo determined by 



tanh(2(/5o) 



Vj — ujv — 2U\/h 
vj + v N + 2U\/h' 



(8) 



Thus, the Luttinger liquid relation survives with the ef- 
fective Luttinger liquid parameter 



A o// =e- 2 *>. 



(9) 



The non-ideal gas theory does not work for the pure 
Lieb-linger gas because eq.(J3J is back to the free bo- 
son when A — * 0. For A = 0, the model is exactly 
soluble. We can use the excluson gas theory with mu- 
tual statistics which was also developed in ref. [l^. It 
is a bosonization theory with a similar Hamiltonian to 
(|7|) with the effective Luttinger liquid parameter X e ff = 
[z(k^)z(kp ^J -1 where z(k) is given by the integration 
equation z{k) = 1 + f% % {Ul/2 ^ k , y , z(k>). The 
pseudo-Fermi momentum kp is determined by the van- 
ishing point of the dressed energy 0,0] and kp — kp±e 



for 



H 



In the T-G limit {U\ — > oo), one hnds 



z(k { p ] ) = 1 and z(k F ] ) = 2, i.e., A e // = 1/2 fl. 

Now, we discuss the possible experimental implica- 
tion. Xeff — 1 is the exponent of the one-dimensional 
momentum distribution, which can be measured by the 
absorption image of the time-flight. Keeping the inner 
wall potential of the cylinder and switching off the outer 
wall potential, the atom cloud will freely expand. The 
atoms with the fastest radial velocity will run along the 
blue lines in Fig. 03 if they have the different azimuthal 
momenta p. Counting the number of the different arc 
intervals (pt) at a given time t, one can draw out the mo- 
mentum distribution of the cigar-like condensates before 
releasing. Then the exponent A e // can be determined. 
The system with A ^ is very different from A = 0. 
This can be seen from the limit U\ — ► oo. For the for- 



There is a restriction in for a given A because ^ A ^ ^ and the system is an insulator while the 



|tanh20 o | < 1. For A > 0, this holds if U\ > -\Xhv^ 
Fig. ^ shows the Luttinger liquid parameter changes as 
interaction. An attractive U\ suppresses A e // from A 
while a repulsive one lifts A e // as expected. This ex- 
ponent X e f f may determine the non- universal exponents 
in the Luttinger liquid |22|. For example, the asymp- 
totic single particle correlation function is proportional 



eff 

latter has A 



eff 



1/2, the T-G gas state. 
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